skip to main content


Search for: All records

Creators/Authors contains: "Liu, Michael C."

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. Abstract

    We present a uniform forward-modeling analysis of 90 late-M and L dwarfs in nearby young (∼10–200 Myr) moving groups, the Pleiades, and the Hyades using low-resolution (R≈ 150) near-infrared (0.9–2.4μm) spectra and the BT-Settl model atmospheres. We derive the objects’ effective temperatures, surface gravities, radii, and masses by comparing our spectra to the models using a Bayesian framework with nested sampling and calculate the same parameters using evolutionary models. Assuming the evolutionary-based parameters are more robust, our spectroscopically inferred parameters from BT-Settl exhibit two types of systematic behavior for objects near the M-L spectral type boundary. Several objects are clustered aroundTeff≈ 1800 K andlogg5.5dex, implying impossibly large masses (150–1400MJup), while others are clustered aroundTeff≳ 3000 K andlogg3.0dex, implying unphysically low masses and unreasonably young ages. We find the fitted BT-Settl model spectra tend to overpredict the peakJ- andH-band flux for objects located near the M-L boundary, suggesting the dust content included in the model atmospheres is insufficient to match the observations. By adding an interstellar medium–like reddening law to the BT-Settl model spectra, we find the fits between models and observed spectra are greatly improved, with the largest reddening coefficients occurring at the M-L transition. This work delivers a systematic examination of the BT-Settl model atmospheres and constitutes the largest spectral analysis of benchmark late-M- and L-type brown dwarfs to date.

     
    more » « less
  2. Abstract

    We derive the bolometric luminosities (Lbol) of 865 field-age and 189 young ultracool dwarfs (spectral types M6–T9, including 40 new discoveries presented here) by directly integrating flux-calibrated optical to mid-infrared (MIR) spectral energy distributions (SEDs). The SEDs consist of low-resolution (R∼ 150) near-infrared (NIR; 0.8–2.5μm) spectra (including new spectra for 97 objects), optical photometry from the Pan-STARRS1 survey, and MIR photometry from the CatWISE2020 survey and Spitzer/IRAC. OurLbolcalculations benefit from recent advances in parallaxes from Gaia, Spitzer, and UKIRT, as well as new parallaxes for 19 objects from CFHT and Pan-STARRS1 presented here. Coupling ourLbolmeasurements with a new uniform age analysis for all objects, we estimate substellar masses, radii, surface gravities, and effective temperatures (Teff) using evolutionary models. We construct empirical relationships forLbolandTeffas functions of spectral type and absolute magnitude, determine bolometric corrections in optical and infrared bandpasses, and study the correlation between evolutionary model-derived surface gravities and NIR gravity classes. Our sample enables a detailed characterization ofBT-SettlandATMO2020 atmospheric model systematics as a function of spectral type and position in the NIR color–magnitude diagram. We find the greatest discrepancies between atmospheric and evolutionary model-derivedTeff(up to 800 K) and radii (up to 2.0RJup) at the M/L spectral type transition boundary. With 1054 objects, this work constitutes the largest sample to date of ultracool dwarfs with determinations of their fundamental parameters.

     
    more » « less
  3. Abstract

    We analyzed 20 s cadence Transiting Exoplanet Survey Satellite time-series photometry of the exoplanet host star HR 8799 collected in Sector 56. The amplitude spectrum shows Gamma Doradus (γ Dor) pulsations consistent with previous space-based photometry from MOST. Assuming that HR 8799 is a representative ofγ Dor stars in the Kepler sample, the dominant dipole mode at 1.98 cycles day−1implies a core rotation period of ∼0.7 day, which combined withvsiniand stellar radius measurements would result in a preliminary stellar inclination of ∼28° assuming rigid rotation. We find no significant residual photometric variation after removing the pulsation signal aside from a ∼9 days trend that is likely a systematic effect or an artifact from performing aggressive frequency subtraction in the presence of red noise.

     
    more » « less
  4. Abstract

    51 Eri is well known for hosting a directly imaged giant planet and for its membership to theβPictoris moving group. Using 2 minute cadence photometry from the Transiting Exoplanet Survey Satellite (TESS), we detect multiperiodic variability in 51 Eri that is consistent with pulsations of Gamma Doradus (γDor) stars. We identify the most significant pulsation modes (with frequencies between ∼0.5 and 3.9 cycles day−1and amplitudes ranging between ∼1 and 2 mmag) as dipole and quadrupole gravity modes, as well as Rossby modes, as previously observed in KeplerγDor stars. Our results demonstrate that previously reported variability attributed to stellar rotation is instead likely due toγDor pulsations. Using the mean frequency of the= 1 gravity modes, together with empirical trends of the KeplerγDor population, we estimate a plausible stellar core rotation period of0.90.1+0.3days for 51 Eri. We find no significant evidence for transiting companions around 51 Eri in the residual light curve. The detection ofγDor pulsations presented here, together with follow-up observations and modeling, may enable the determination of an asteroseismic age for this benchmark system. Future TESS observations would allow a constraint on the stellar core rotation rate, which in turn traces the surface rotation rate, and thus would help clarify whether or not the stellar equatorial plane and orbit of 51 Eri b are coplanar.

     
    more » « less
  5. Abstract We present the third discovery from the COol Companions ON Ultrawide orbiTS (COCONUTS) program, the COCONUTS-3 system, composed of the young M5 primary star UCAC4 374−046899 and the very red L6 dwarf WISEA J081322.19−152203.2. These two objects have a projected separation of 61 ′ ′ (1891 au) and are physically associated given their common proper motions and estimated distances. The primary star, COCONUTS-3A, has a mass of 0.123 ± 0.006 M ⊙ , and we estimate its age as 100 Myr to 1 Gyr based on its stellar activity (via H α and X-ray emission), kinematics, and spectrophotometric properties. We derive its bulk metallicity as 0.21 ± 0.07 dex using empirical calibrations established by older and higher-gravity M dwarfs and find that this [Fe/H] could be slightly underestimated according to PHOENIX models given COCONUTS-3A’s younger age. The companion, COCONUTS-3B, has a near-infrared spectral type of L6 ± 1 int-g , and we infer physical properties of T eff = 1362 − 73 + 48 K, log ( g ) = 4.96 − 0.34 + 0.15 dex, R = 1.03 − 0.06 + 0.12 R Jup , and M = 39 − 18 + 11 M Jup using its bolometric luminosity, its host star’s age, and hot-start evolution models. We construct cloudy atmospheric model spectra at the evolution-based physical parameters and compare them to COCONUTS-3B’s spectrophotometry. We find that this companion possesses ample condensate clouds in its photosphere ( f sed = 1) with the data–model discrepancies likely due to the models using an older version of the opacity database. Compared to field-age L6 dwarfs, COCONUTS-3B has fainter absolute magnitudes and a 120 K cooler T eff . Also, the J − K color of this companion is among the reddest for ultracool benchmarks with ages older than a few hundred megayears. COCONUTS-3 likely formed in the same fashion as stellar binaries given the companion-to-host mass ratio of 0.3 and represents a valuable benchmark to quantify the systematics of substellar model atmospheres. 
    more » « less
  6. ABSTRACT

    VHS J1256−1257 AB is an ultracool dwarf binary that hosts a wide-separation planetary-mass companion that is a key target of the JWST Exoplanet Early Release Science programme. Using Keck adaptive optics imaging and aperture masking interferometry, we have determined the host binary’s orbit (a = 1.96 ± 0.03 au, P = 7.31 ± 0.02 yr, e = 0.883 ± 0.003) and measured its dynamical total mass (0.141 ± 0.008 M⊙). This total mass is consistent with VHS J1256−1257 AB being a brown dwarf binary or pair of very low-mass stars. In addition, we measured the orbital motion of VHS J1256−1257 b with respect to the barycentre of VHS J1256−1257 AB, finding that the wide companion’s orbit is also eccentric ($e=0.68^{+0.11}_{-0.10}$), with a mutual inclination of 115○ ± 14○ with respect to the central binary. This orbital architecture is consistent with VHS J1256−1257 b attaining a significant mutual inclination through dynamical scattering and thereafter driving Kozai–Lidov cycles to pump the eccentricity of VHS J1256−1257 AB. We derive a cooling age of 140 ± 20 Myr for VHS J1256−1257 AB from low-mass stellar/substellar evolutionary models. At this age, the luminosity of VHS J1256−1257 b is consistent with both deuterium-inert and deuterium-fusing evolutionary tracks. We thus find a bimodal probability distribution for the mass of VHS J1256−1257 b, either 12.0 ± 0.1 MJup or 16 ± 1 MJup, from these models. Future spectroscopic data to measure isotopologues such as HDO and CH3D could break this degeneracy and provide a strong test of substellar models at the deuterium-fusion mass boundary.

     
    more » « less
  7. Abstract

    M dwarfs are common host stars to exoplanets but often lack atmospheric abundance measurements. Late-M dwarfs are also good analogs to the youngest substellar companions, which share similarTeff∼ 2300–2800 K. We present atmospheric analyses for the M7.5 companion HIP 55507 B and its K6V primary star with Keck/KPIC high-resolution (R∼ 35,000)K-band spectroscopy. First, by including KPIC relative radial velocities between the primary and secondary in the orbit fit, we improve the dynamical mass precision by 60% and findMB=88.03.2+3.4MJup, putting HIP 55507 B above the stellar–substellar boundary. We also find that HIP 55507 B orbits its K6V primary star witha=383+4au ande= 0.40 ± 0.04. From atmospheric retrievals of HIP 55507 B, we measure [C/H] = 0.24 ± 0.13, [O/H] = 0.15 ± 0.13, and C/O = 0.67 ± 0.04. Moreover, we strongly detect13CO (7.8σsignificance) and tentatively detectH218O(3.7σsignificance) in the companion’s atmosphere and measure12CO/13CO=9822+28andH216O/H218O=24080+145after accounting for systematic errors. From a simplified retrieval analysis of HIP 55507 A, we measure12CO/13CO=7916+21andC16O/C18O=28870+125for the primary star. These results demonstrate that HIP 55507 A and B have consistent12C/13C and16O/18O to the <1σlevel, as expected for a chemically homogeneous binary system. Given the similar flux ratios and separations between HIP 55507 AB and systems with young substellar companions, our results open the door to systematically measuring13CO andH218Oabundances in the atmospheres of substellar or even planetary-mass companions with similar spectral types.

     
    more » « less
  8. Abstract

    We present the characterization of the low-gravity M6 dwarf 2MASS J06195260-2903592, previously identified as an unusual field object based on its strong IR excess and variable near-IR spectrum. Multiple epochs of low-resolution (R≈ 150) near-IR spectra show large-amplitude (≈0.1–0.5 mag) continuum variations on timescales of days to 12 yr, unlike the small-amplitude variability typical for field ultracool dwarfs. The variations between epochs are well-modeled as changes in the relative extinction (ΔAV≈ 2 mag). Similarly, Panoramic Survey Telescope and Rapid Response System 1 optical photometry varies on timescales as long as 11 yr (and possibly as short as an hour) and implies comparableAVchanges. Near Earth Object Wide-field Infrared Survey Explorer mid-IR light curves also suggest changes on 6 month timescales, with amplitudes consistent with the optical/near-IR extinction variations. However, near-IR spectra, near-IR photometry, and optical photometry obtained in the past year indicate that the source can also be stable on hourly and monthly timescales. From comparison to objects of similar spectral type, the total extinction of 2MASS J0619-2903 seems to beAV≈ 4–6 mag, with perhaps epochs of lower extinction. Gaia Early Data Release 3 (EDR3) finds that 2MASS J0619-2903 has a wide-separation (1.′2 = 10,450 au) stellar companion, with an isochronal age of3110+22Myr and a mass of0.300.03+0.04M. Adopting this companion’s age and EDR3 distance (145.2 ± 0.6 pc), we estimate a mass of 0.11–0.17Mfor 2MASS J0619-2903. Altogether, 2MASS J0619-2903 appears to possess an unusually long-lived primordial circumstellar disk, perhaps making it a more obscured analog to the “Peter Pan” disks found around a few M dwarfs in nearby young moving groups.

     
    more » « less
  9. Abstract

    In order to assess the multiplicity statistics of stars across spectral types and populations in a volume-limited sample, we censused nearby stars for companions with Robo-AO. We report on observations of 1157 stars of all spectral types within 25 pc with decl. >−13° searching for tight companions. We detected 154 companion candidates with separations ranging from ∼0.″15 to 4.″0 and magnitude differences up to Δmi7using the robotic adaptive optics instrument Robo-AO. We confirmed physical association from Gaia EDR3 astrometry for 53 of the companion candidates, 99 remain to be confirmed, and two were ruled out as background objects. We complemented the high-resolution imaging companion search with a search for comoving objects with separations out to 10,000 au in Gaia EDR3, which resulted in an additional 147 companions registered. Of the 301 total companions reported in this study, 49 of them are new discoveries. Out of the 191 stars with significant acceleration measurements in the Hipparcos–Gaia catalog of accelerations, we detect companions around 115 of them, with the significance of the acceleration increasing as the companion separation decreases. From this survey, we report the following multiplicity fractions (compared to literature values): 40.9% ± 3.0% (44%) for FGK stars and 28.2% ± 2.3% (27%) for M stars, as well as higher-order fractions of 5.5% ± 1.1% (11%) and 3.9% ± 0.9% (5%) for FGK stars and M-type stars, respectively.

     
    more » « less